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ABSTRACT 

I pinpoint a previously unrecognized MOND effect that may act in the inner solar 
system, and is due to the galactic acceleration, gg = r/Oo: a byproduct of the MOND 
external- field effect. Predictions of the effect are not generic to the MOND paradigm, 
but depend on the particular MOND formulation at hand. However, the modified- 
Poisson formulation, on which I concentrate, uniquely predicts a subtle anomaly that 
may be detected in planetary and spacecraft motions (and perhaps in other precision 
systems, such as binary pulsars), despite their very high accelerations, and even if the 
MOND interpolating function is arbitrarily close to unity at high accelerations. Near 
the sun, this anomaly appears as a quadrupole field, with the acceleration at position 
u from the sun being gf(u) = —q^jU^, with q'^j diagonal, axisymmetric, and traceless: 
-2q^^ = -2q^y = q^^ = q{ri){a^ j R^) , where Rm = (M^G/a^y^ « 8 • lO^a.u. is the 
MOND radius of the sun. The anomaly is described and analyzed as the Newtonian 
field of the fictitious cloud of "phantom matter" that hovers around the sun. I find, for 
the relevant range of ry values, and for a range of interpolating functions, values 
of 10^^ < —q < 0.3, which turn out to be sensitive to the form of ^{x) around the 
MOND-to-Newtonian transition. This range verges on the present bounds from solar 
system measurements. There might thus exist favorable prospects for either measuring 
the effect, or constraining the theory and the relevant parameters. Probing this anomaly 
may also help distinguish between modified inertia and modified gravity formulations of 
MOND. I also discuss briefly an anomaly that is generic to MOND in all its formulations, 
and competes with the quadrupole anomaly in the special case that 1 — n{x) vanishes 
as as X — > oo. 

Subject headings: 



1. Introduction the MOND external- field effect 

The MOND external-field effect (EFE) encapsulates the fact that even a constant-acceleration 
gravitational field in which a system is falling may enter and affect the internal dynamics of this 
system. It results from the non-linearity of MOND, which, in turn, follows from the basic premises 
of the paradigm. This means that even though all the constituents of the system fall in the same 
external field, this field does not cancel in their relative motion, as happens in (the linear) Newtonian 
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dynamics. The effect appears already in the pristine formulation of MOND (Milgrom 1983), as 
well as in the formulation of MOND as a modified gravity described by a nonlinear extension of the 
Poisson equation (Bekenstein & Milgrom 1984, Milgrom 1986a). Many applications to astrophysical 
systems have been discussed to date, including dwarf spheroidal galaxies in the field of a mother 
galaxy (Brada k, Milgrom 2000a, Angus 2008), globular clusters (Haghi et al. 2009), warp induction 
by a companion (Brada k, Milgrom 2000b), escape speed from a galaxy (Famaey, Bruneton, k, Zhao 
2007, Wu et al. 2007), binary galaxies (Tiret et al. 2007), departure from asymptotic flatness of 
the rotation curve (Wu et al. 2007), galaxy clusters (Angus k, McGaugh 2008), and others. 

In all these applications one deals with accelerations, both internal (gj„) and external (gea,), that 
are of order of, or smaller than, the MOND acceleration constant Oq, and with external accelerations 
that are similar or larger than the internal ones. In these cases the pristine formulation and the 
modified-Poisson equation give similar results, barring order-of-unity differences in the geometry. 
In modified-inertia formulations of MOND the EFE still acts, but its implications can be quite 
different (Milgrom, in preparation). 

Here I consider a situation in which the relative magnitudes of the three accelerations are very 
different: gi„ gg > Uq (in what follows I use everywhere gex = Eg for the galactic field). This 
characterizes, for example, high acceleration systems in the galactic neighborhood of the sun, such 

as the inner solar system (SS) or binary stars. These three inequalities would each, separately, 
reduce the effect, and render it small. However, the effect might still be detectable in systems that 
are measured with high accuracy such as motions in the solar system, or in binary pulsars. This 
arrangement of accelerations does not lend itself to simple analytic, or numerical, considerations, 
as do the cases studied so far. 

Small, nonrelativistic MOND anomalies in high-acceleration motions may possibly have several 
origins. In non-local, modified inertia formulations of MOND, anomalies may appear in the inner 
solar system in the motion of bodies that are on highly eccentric trajectories; trajectories that take 
them to large distances, where accelerations are low (Milgrom, in preparation). Cases in point are 
the long-period comets, and the Pioneer spacecraft. Because of the non-locality of such theories, 
the motions in high-acceleration portions of the orbit may be affected by the low acceleration ones, 
giving rise to small anomalies. Such MOND effects have been proposed as a possible mechanism for 
generating the Pioneer anomaly, without affecting the motions of planets, whose orbits are wholly 
in the high acceleration regime (Milgrom 2002, 2005, in preparation). 

I concentrate on nonrelativistic, modified-gravity formulations, in which MOND is described 
by a modified potential field, in which particles move according to the standard laws of motion. 
In this case, an anomaly can arise in high acceleration systems due to a remaining departure of 
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the MOND interpolating function, n{x), from unity, even at high accelerations.^'^ For example, in 

the inner SS such a departure appears as an anomalous, spherical acceleration field. Its effects on 
planetary motions, and constraints on the high-x behavior of 5fi{x) = 1 — fi{x), have been discussed, 
e.g., in Milgrom (1983), and in Sereno &: Jetzer (2006). In such analyses the sun is treated as an 
isolated point mass. 

Here, I discuss another type of anomaly having to do with possible influence of the galactic 
field through the EFE. This, as we shall see, appears in the inner SS as an anomalous quadrupole 
field, not as a spherical one. 

In order to isolate this effect from the first type, I shall use, all along, forms of with Sfi^x) 
fast enough to dwarf the first effect in comparison with the EFE. 

It is worth noting first that for the parameter values relevant to the SS, there is no unique 
MOND prediction of the EFE. The primitive, pristine formulation (Milgrom 1983), in which the 
MOND acceleration, g, is calculated from the local Newtonian value, g^^, via an algebraic relation 
of the form 

g = i^(|gjv|/ao)giv (1) 

[where yfiy) = x is the inverse of xii{x) = y], captures the main MOND effects in many instances, 
for example, in describing galaxy rotation curves, or the EFE when gg ^> (Milgrom 1983). 
However, it is quite inapplicable to the case Qg ^ such as in solar system, as it is known not 
to give correctly the motion of such compound systems in an external field: For example, atoms in 
a star, or the stars in a close binary, moving in the outskirts of a galaxy are subject to high total 
accelerations because of the internal forces in the system. The naive formulation will say then that 
the full motion of these constituents is Newtonian. But this would imply that the whole star (or 
binary) is also moving Newtonically in the galaxy, contrary to what is observed, and to what we 
want MOND to accomplish. This is also the situation for the solar system, as the acceleration of 
the Sun due to Jupiter alone is ~ 2500ao, and that of all the planets due to the sun, even larger. 
Applying the algebraic formulation to the motions of the Sun and planets, separately, would give 
not only the wrong motion for the solar system in the Galaxy, but also a negligible EFE for a fast 
decreasing dfi{x). 

Modified-inertia formulations of MOND are more versatile and it is not possible to deduce 
a generic modified inertia prediction of the EFE (Milgrom, in preparation). For example, such 
formulations can be constructed in which there is no EFE in the inner SS. This is the case, e.g., in the 



^Galactic dynamics probe this function only at low values of a; < 5, and gives no constraint whatever on the high 
X behavior. It is easy to invent forms of with slow rise near a; = 1, but very fast convergence to 1 at high x (see 
some examples in Milgrom & Sanders 2008). 

^Anomalies may also appear in relativistic theories such as TcVcS and the like, because these theories do not 
exactly coincide with general relativity even in the limit Oo — > (see Bekenstein 2006, Sanders 2006, 2007, Sagi 2009, 
and Skordis 2009.) 
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toy theory described in Milgrom (1999). (This theory does exhibit an EFE in other circumstances.) 
But we do not yet have an acceptable formulation of MOND as modified inertia. 

However, it turns out that in the modified Poisson formulation of Bekenstein &; Milgrom (1984) 
there is an EFE that is felt in the inner SS as a quadrupolc anomaly (to dominant order). This 
anomaly exists despite the high accelerations of the planets, and, interestingly, it remains finite 
even with 6ijl{x) arbitrarily fast at high x. Furthermore, the strength of the effect is sensitive 
to the form of ijl{x) in the MOND-to- Newtonian transition region a; ~ 1, even though the effect 
itself is probed at very high accelerations.^ 

All this results from properties of the Poisson equation, and its ilk, such as the non-linear 
version treated here: Distortion of the field in one place is felt in all other regions. Here, the 
distortion is created in the sub Oq region that must always exist in the direction of the galactic 
center, due to the tag-of-war between the Sun and the Galaxy. 

My main purposes here are to point out the anomaly, describe its nature, calculate its strength, 
and describe some of its effects on planetary motions. I leave more detailed and systematic study 
to future work. 

In section 2, I describe the anomaly in the frame of the modified Poisson theory. In section 
3, 1 briefly comment on multi-potential theories, such as the nonrelativistic limit of TcVcS. In 
sections 4, 1 describe an approximation for computing the effect. In section 5, I give the numerical 
results for various forms of //(x) both for numerical solutions of the full field equation and for the 
approximate expression. In section 6, I look at solar system effects, mainly the induced anomalous 
perihelion precession. Section 7 is a discussion. The appendices describe the properties of the 
fictitious "phantom matter" that might be viewed as giving rise to the anomaly within a Newtonian 
framework. 

2. The quadrupole anomaly in the modified Poisson theory 

The special case of the EFE studied here affects all high-acceleration galactic systems and 
could, in principle, be detected in precision systems such as the SS, binary pulsars, etc.. I shall 
concentrate here on the EFE for a point mass, because it is simpler to treat, involving, as it 
does, only MG, Oq, and r] = ggja^ as defining parameters. While the principles are similar, the 
treatment of a general few-body system is more complicated, because it involves several additional 
dimensionless parameters, such as characteristic lengths in units of the MOND radius of the system, 
Rm = (MG/oo)^/^, the mass ratios, orientation with respect to the external field, etc.. Even the 



^In Milgrom (1986a) I discussed the anomalous potential difference that is predicted between the sun's vicinity 
and infinity. This, however, does not affect motions near the sun, and depends on system parameters, and on the 
choice of n{x), very differently from the quadrupole anomaly discussed here. For example, for the isolated case, rj = 0, 
the quadrupole vanishes, while the potential anomaly becomes infinite. 
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isolated, two-body problem is not exactly solvable in MOND, while the test particle motion in the 
field of a point mass is. The EFE for few-body systems will be discussed elsewhere. 

Accordingly, I consider here the motion of a test particle (e.g., a planet) in the field of a 
point mass (e.g., the sun), embedded in an asymptotically constant gravitational field. The latter 
approximates, e.g., the galactic field at the position of the sun.^ In this case, the acceleration of 
the test particle is given by the gradient of the modified potential field (not so for a non-test mass); 
so, we set to calculate the anomalous potential field for this problem. 

In the modified Poisson theory of Bekenstein & Milgrom (1984), the gravitational potential, 
4>{r), due to a density p(r), embedded (falling) in a constant external field g^, is gotten from the 
nonlinear generalization of the Poisson equation: 

V • [/x(|V(/>|/ao)V<^] = AirGp. (2) 

Here /x(a:) is the interpolating function characterizing the theory. This equation has to be solved 

with the boundary conditions: V(/> ~Sg- The internal dynamics is then governed by the field 
V(j)i„ = V(/) + gg. The MOND anomaly is that part of Vcpin in excess of the Newtonian field for 
the same mass when isolated, call it — V^^. The MOND anomaly is then 

g" = -V<^a = -(V.^.n-Vr)- (3) 

We want to determine the behavior of this anomaly at distances from the Sun much smaller than 
Rm, where the accelerations are much higher than Cq. 

The above anomalous acceleration describes all MOND effects combined. When 6iJ,{x) = 
1 — iJ,{x) vanishes at high x as Ax~°', we get an anomaly near the sun, even in the absence of 
an external field. This isolated-mass acceleration is spherical and pointing towards the sun, g" = 
—Aao{u/RM)'^^'^~^^e, where e = u/u is the radial unit vector and u is the position with respect 
to the Sun (Milgrom 1983). For example, for a = 1 wc get a constant anomalous acceleration 
pointing towards the sun g" = —Auoe. This could explain the Pioneer anomaly, but produces 
too strong effects on the planets (as already shown by the analysis of Milgrom 1983). In a recent 
detailed analysis, Sereno &; Jetzer (2006) conclude that SS measurements roughly allow only a > 1.5 
(depending on the value of A). 

Sereno & Jetzer take the parameter values (such as perihelion precession rate, or semi-major 
axis) quoted as yielding the minimum best ephemerides fit residuals, as actual limits on any possible 
anomalous force, not included in the fits. They thus obtain limits on the strength of such various 
forces, including MOND anomalies for various a. However, one may question this procedure, since 
the more sensible approach would be to include any suspected anomalous force in a global fit. 



''in principle, we have to solve for the combined Galaxy-sun field. Our approximation amounts to taking the 
Galaxy as a mass M — » oo, at a distance ii — > oo, with M/ii^ fixed. Corrections to this approximation are expected 
to be of order (Mo/M)^/^ ~ 3 • 10"^ 
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and in this way to determine what strength of the anomaly can be tolerated by the data (I am 
grateful to Gillcs Esposito-Farese for making this point to me). For example, using the above 
criterion, Sercno and Jetzer deduce an upper limit of lO^^'^cm s^^ on a constant, anomalous, 
sunward acceleration that can be tolerated by the data for Uranus, 4 • 10"^'' cm s~^ for Neptune, 
and 10~^cm s"^ for Pluto. However, Fienga et al. (2009), following the more sensible procedure, 
and including a constant acceleration in the equations of motion of these three planets in their 
global fitting procedure (to test the tolerance to a Pioneer-like anomaly) , find that the Uranus data 
can tolerate a constant acceleration of up to ~ 2 • 10~^cm s~^, and that of Neptune and Pluto 
allow an anomalous constant acceleration even of the size of the claimed Pioneer anomaly; namely 
~ 8 • 10~^cm s~^. These are some two order of magnitude larger than the limits deduced by Sereno 
and Jetzer (2006). Clearly, the whole question of SS constraints on the high acceleration behavior 
of fx{x) requires revisitation. 

Be that as it may, the case a = 3/2 happens to be interesting conceptually: It gives rise to 
a harmonic, anomalous acceleration g" = —Aaou/RM, which has the same u dependence, and a 
similar magnitude, as the quadrupole anomaly that is my main subject in this paper. 

In most of what what follows I isolate the EFE from such spherical anomalies by assuming 
forms of fi with fast decreasing d^{x) at large values of x. I will, however, discuss, in section 6, 
the effects of the harmonic anomaly produced in the case of the critical exponent a = 3/2, which 
is interesting in the present context. It is an effect that is generic to MOND and obtain in any 
formulation as long as the characteristic interpolation function has the appropriate asymptotic 
form. The two effects are separate and largely additive, since the quadrupole anomaly depends 
mainly on the behavior of fi near the transition region, while the harmonic anomaly derives from 
the high x behavior of fx{x). 

What then is g" like near the Sun under the above restriction on fi7 Note first, importantly, 
that g" vanishes identically at the position of the Sun (shown in Appendix A.l). Secondly, note that 
with our choice of fJ,{x), the anomaly is practically harmonic (a solution of the Laplace equation) 
within a large volume around the sun. The reason for this is as follows: in this region |V(/)|/a(, ^ 1 
and so, practically, jj, = 1, there. As a result, (p solves the (linear) Poisson equation there, with 
the same density source (the sun) as 0^, and hence their difference is harmonic. (For a < 3/2 
this is not the case: for Example, for a = 3/2 we have a constant phantom density near the sun.) 
Add to these the fact that the only length scale in the problem is Rm ~ 8 • lO^a.u., and that the 
dimensionless parameter 77 is of order unity, and we conclude that when |u| is small enough with 
respect to Rm, '^a(u) is dominated by a quadrupole (a dipole is absent because the anomalous 
acceleration vanishes at u = 0); so we can write there 



where the anomalous quadrupole is diagonal, and axisymmetric in a frame where the z axis is along 
the external field. Harmonicity of (f)a implies that g"- is traceless: 



</>a(u) «i .^a(O) + -qfjU'u^, 



(4) 




= -9^/2. 



(5) 
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The effect is tlius cliaracterized by one quantity tfiat we fiave to determine, say q^^. On dimensional 
grounds we can write: 

= q{v)a,R-\ (6) 

and I seek to determine tlie dimensionless, quadrupole parameter q{r]), wfiicfi describes tfie EFE 
anomaly.^ 

Near the sun, the anomaly, as defined in eq.(3), is a very small difference between the much 
larger V0„, and V^'^. We have |V(/)a|/| Vc/)'^] = \q\{u/ RM)[aa/ {MG /u^)]\ the last two factors are 
very small (their product is ~ 10~^^ at la.u. from the sun), and we shall see that q ~ 0.1. Deducing 
the anomaly directly from eq.(3) by solving numerically for cj) (e.g., using a non-linear Poison solver) 
could be very challenging. Now that we identified the position dependence of the anomaly, one can 
perhaps deduce q by evaluating the anomaly numerically at a larger distance from the Sun (e.g., 
at lO^a.u. the relative strength of the anomaly is ~ 10~^). 

Even though the effect is very small in the inner solar system, the problem cannot be attacked 

by expanding the governing equations in a small perturbation: There is always a region where the 
MOND field differs greatly from the Newtonian one (sec below), and the ellipticity of the modified 
Poisson equation ties the different regions inseparably. 

There is a useful way to envisage the anomaly based on our intuition from Newtonian gravity, 
using the concept of "phantom mass" introduced in Milgrom (1986b). This is the density distri- 
bution, Pp, whose Newtonian field is the MOND anomaly. It is the mass density that a Newtonist 
would attribute to dark matter, and is given by 

" ^^^^ = • - = 4^^ • - 

We can also write 

= • W - p = • [(1 - /x)V<^]. (8) 

The exact MOND field in the system, —Vcp, equals the sum of the Newtonian fields of p and pp, 
and of the constant (MOND) external field. ^ 

In this way, wc have identified the MOND anomaly as the full field of a certain mass dis- 
tribution, instead of as a very small difference between two strong fields. This can be useful for 
deducing the anomaly numerically, and for devising approximation schemes that will bypass the 
inherent difficulty of deducing a small difference of large quantities. 



^Since the EFE anomalous acceleration is linear in u, it dominates the first type of anomaly, near the sun, when 
S/j, oc x"" at high x, for a > 3/2. 

®Note that this field differs from the field of a central quadrupole, such as an oblate sun, which has the same 
angular dependence but a different radial one. 

'^The MOND external field can also be viewed as the sum of the Newtonian galactic field, plus that of the galactic 
phantom density. 
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Using the field equation (2) we can write 



1 



(m7m)v|V0|-w + U 



-1 



(9) 



Pp 



or in another useful form 




We see from expression (9) that pp resides where jU differs from 1, namely, in low acceleration 
regions. There must exist such a region, even when the external field itself is high: Near the Sun the 
acceleration is pointing towards the sun. As we move away in the direction of the galactic center, 
the acceleration must eventually vanish at some point, and change sign to point in the opposite 
direction. There is thus always a region V around the critical point that is in the sub a„ regime, 
in which MOND acts in full force, and in which pp is not negligible. This region excludes a large 
volume around the Sun with our choice of p,. Even with the limiting form of p, where = 1 for 
X > 1, the region V, the distribution of pp in it, and thus the anomaly that is felt near the sun, 
remain finite. 

The anomaly can be written as 



where u = r — Rq is the position from the sun, and cosO = Uz = Uz/u. 

In Appendix A, I derive several useful, exact properties of the phantom mass. First, I show that 
around the critical point there are always regions of both positive and negative phantom density 
(A.3) (this is an example of the general trend discussed in Milgrom 1986b). Second, in A.4 I derive 

the total phantom mass Mp; it depends on rj and the form of p. (For the limiting form of p we 
have Mp = identically; so the negative and positive densities exactly cancel) . Third, I derive the 
column density along the symmetry axis (A. 2). It is independent of ry, and for the limiting form of 
p it is ao/27rG. Also, I find [in A. 6] that pp decreases asymptotically as (taking both signs at 
all radii); so this provides good convergence in eq.(12). 



This can be used only when pp vanishes at the origin, as I assume all along. Otherwise the u integration diverges 
at the origin and the integral is undefined. The quantity this integral describes can still be finite when pp doesn't 
vanish at the origin. 




(11) 



We thus identify^ 




(12) 
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I also derive in the appendix [eq.(A12)] an expression for the anomaly as an integral over 
(not involving derivatives of V(f) as appear in pp): 

where /x^ = //(r;) is the galactic value of jj., and n is a unit vector in the direction of r — R. Although 
eq.(13) is equivalent to eq.(3), it is a better starting point for approximations. The great sensitivity 
of the anomalous acceleration, as deduced from cq.(3), to errors in (or equivalently in V(^i„), 
still remain in expression (13), but here it is controllable. Errors in the integration, will cause 
g"(R) to vanish not at the position of the sun, but at a somewhat different point off by (5R. Even 
if |(5R|/i?M ^ 1) the resulting errors in g"(R) near the Sun could be very large. However, we can 
eliminate this problem by considering directly^ 



dZ 



= - J- / ^l-u^ N.-^- + 3V(/) • n(l - 5nl)]d\ (14) 
. 47r 7 |u|^ 



where n is now a unit vector in the direction of u. If in eq.(14) we take u in units of Rm and V0 
in units of Oo, the expression yields q. 

The phantom mass is only a visualization aid that I find very convenient. We could have 
derived eqs.(13) and (14), simply by starting from the identity (for functions that decay fast enough 
at infinity): 

(where the second equality comes from the definition of 0^), and proceed from there, integrating by 
parts, without mentioning phantom matter. The concept is, however, also useful because my main 
approximation is to replace the true phantom matter by a surrogate one, which I do in section 4. 

To recapitulate, eqs.(ll)-(14) are only identities, useful here in two ways: a. for calculating the 
quadrTipole i^"^ as an integral over "large" quantities, instead of as a very small difference between 
much larger quantities, b. for devising approximations for g^^. 

To use eq.(14) to calculate the anomaly we need to solve the nonlinear Poisson equation for 
the required choice of interpolating function and an rj value, and use the resulting field in eq.(14) 
(instead of subtracting the Newtonian acceleration from the numerical solution very near the sun) . 



®This expression converges well near the Sun despite the |u|* factor, and the divergence of V0, because A* — 1 
vanishes quickly there (and because the field becomes spherical, for which case the angular integral vanishes). When 
Sfi{x) vanishes faster than x~^^^ for x ^ oo, the u integration is convergent near « = 0, and we can use this 
expression without further care. For 6^{x) x~^^^, the u integral diverges logarithmically; however, the angular 
integral vanishes, and the z derivative of the force is still finite at the sun. In this case, the anomalous field near 
the sun is not the pure quadrupole, as I mentioned above. At infinity, the convergence is also good, because there 
— > — gg + V(5 with V5 decreasing like u^^. For the constant part the angular integral vanishes, so there is 
convergence at least as d^u/u^. See also eq.(12) 
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It is useful, however, to develop an approximation scheme that will permit a quick calculation of q 

for mass-production studies, since /j, and the relevant i] for the solar system are not exactly known, 
and also for application to other precision systems, perhaps involving more massive bodies. I shall 
then resort to an approximation that will replace pp with a reasonable surrogate, pp, that is known 
in closed form. The surrogate anomaly can then be written as an integral of a known function. 



3. Multi-potential theories 

It is interesting to study the EFE in the SS in other modified- gravity formulations of MOND, 
such as multi-potential, modified-Poisson theories. In these, the gravitational potential is the sum 
of several potentials, each satisfying a modified Poisson equation of its own. For example, the 
nonrelativistic limit of TeVeS is a double-potential theory: the potential is the sum of a Poissonian 
potential and one that is a solution of the nonlinear Poisson equation, both having the true density 
as source. Generalizing this, consider a nonrelativistic theory in which a particle's equation of 
motion is a = — V$, with $ = (jji, and with determined from 

V • [pi(\V(ki\/aia,)Vct)i] = 47rAGp. (16) 

For concreteness sake, assume that in the limit x ^ oo all pi{x) go to constants^'^ , and take (3i 
such that all these constants are 1. (See relevant comments in Zhao & Famaey 2006, Famaey & al. 
2007, Bruneton & esposito-Farese 2007, and Milgrom &; Sanders 2008, who deal with constraints 
on the form of the interpolating function in TeVeS.) Also take G to be the usual gravitational 
constant, which means that Ylif^i — 1- The dimensionless coefficients can be absorbed into the 
definition of pi, but it is convenient to keep them so that we can require a uniform \ow-x behavior 
Pi{x) X. In this specific case, brought as an example, we have in the very- low-acceleration limit 
\V4>i\ <C aiOo, for a spherical system, |V0i| ~ (ai/3j)"'^/^(ao| V^^l)"*^/^, and Oq then plays its usual 
role in phenomenology if we normalize it so that = 1 (I assume that ai,l3i > 0, so all 

potentials describe attractive gravity). If one of the CKj — ^ 0, the corresponding pi{x) — > 1 for all 
x; so, (pi satisfies the linear Poisson equation V • V^j = AiifiiGp. For example, the nonrelativistic 
limit of TeVeS would correspond to ai — > 0, and p2 = P and f32 = P still unspecified (/3i = 1 — (3, 
a2 = /?-!). 

Each of the sub-potentials 4>i is solved for separately, each generates its own cloud of "phantom 
matter" (which do not interact with other clouds), each creating its own sub-anomaly 

g« = -v(</.i-An (17) 

(^^ is the standard Newtonian potential, calculated with the standard G). 

To apply our results to the multi-potential case, we first have to decompose the galactic field 
at the system's (e.g., the sun's) position to its different components (for example, by solving the 



In principle we can have more general behaviors, but I do not wish to delve into a general discussion here. 
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full problem for the Galaxy's mass distribution, determining the values of the sub-fields at that 
position): = X^jg^. Then, remembering that the effective constants for each potential are 
Gi = PiG, al = aiGo, and R\j = {(3i/ai)^^'^RM, we define r]i = |g*|/Q;ja(,. We then calculate the 
dimensionless quadrupole anomaly qi{r]i) using our results above and below. We can also, as an 
approximation that is exact for a spherical galaxy, use the algebraic relations 

gp = Ai^(A|g^lMao)g^, (18) 

to determine g* from the Newtonian galactic field g^. 

The Sun is at a position in the Galaxy with a symmetric perspective. Under these simplifying 
circumstances, all the g* must point approximately in the same direction: towards the galactic 
center. The anomaly's symmetry axis is then common to all the sub anomalies, and we simply 
have, again, a biaxial, quadrupole anomaly with — 2g°^ = —2qyy = = qao/RM, and 

q = J2qi{Vi)aT/3;'/'. (19) 
i 

For a planetary system placed at an arbitrary position in the Galaxy, the directions of g* 
may be different. We then straightforwardly add up the different quadrupole sub anomalies. The 
result is still a traceless quadrupole; but it is, in general, triaxial, with all three principal axes, and 
elements, determined by the qi{r]i) and the respective directions of g* . 

For TeVeS, the Poissonian field, (f)i, does not create an anomaly. The galactic field g^ defines 
the symmetry axis, even for an arbitrary location in the Galaxy. The anomaly is thus always a 
biaxial quadrupole, with q = (5~'^q2{V2)- 



4. The surrogate phantom density and its properties 

As an approximation, I shall substitute for pp a surrogate phantom density, pp, derived from 
the divergence of the MOND field, as approximated by the algebraic relation eq.(l). Call this field 
g*, then 

with g* = i^(|giv|/ao)gjv, and 

^qG , u , 

gjv = -riNttoez - -r-nru = -ao{r]j^ez + 7x7^), (21) 
1^1 |u| 

where u = vl/Rm (the galactic field is in the —z direction), and rji^ is the Newtonian value of rj. 
We thus have pp in closed form.^^ 



^^Note that I do not use g* as the MOND field; i.e., I do not substitute it for V0 for use in eq.(3). I only use 
it to obtain an approximation for pp. The Newtonian field of pp; i.e., the surrogate anomaly g", which is derivable 
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All the results based on the surrogate density depend on the Newtonian galactic field at the 
sun's position, g^, and do not require knowledge of the MOND value of this field, g^, while 
the results using the exact theory depend on the latter alone. The relation between these two 
accelerations depends on the galactic mass distribution, which we don't know exactly (we can 
measure gg, but need a mass model to determine g^). We do know that because the Sun is very 
near the galactic mid-plane, symmetry dictates that the two are very nearly parallel. I assume 
in all that follows that they are, and that their values are related by ij^ = m(^)^) which is exact 
for a spherical galaxy. This is a good approximation near the sun, and, in any event, has little 
consequence for us here, since it is only used to relate the value of r/jv to be used in the surrogate- 
density approximation to the value of rj appropriate for the sun: the small uncertainty in rjp, (rj) is 
rater immaterial. 

The surrogate phantom density pp is somewhat different from pp. However, pp does have 
properties (demonstrated in Appendix B) that are rather similar to those of pp, on which we can 
found the hope that this is a reasonable approximation for calculating the anomaly: First, pp has a 
total mass similar to that of pp (for the limiting form of p they both vanish) (B.4). Second, pp has 
a similar pattern of sign distribution (B.3). Third, the column density of pp along the symmetry 
axis is identical to that of pp (B.2). Fourth, the Newtonian field of jjp, like that of pp, vanishes at 
the position of the sun, and is thus also a traceless, quadrupole near the Sun (B.l). 

In some regards, the surrogate density does not quite mimic pp. For example, as can be 
seen from appendices A. 6, and B.6, the angular distributions of their asymptotic behavior (for 
m) are different. 

As shown in A. 5, substituting pp for pp is tantamount to replacing in eq.(13) by — g*. This 
gives for the surrogate anomaly: 

r(R) = ^^[Sn - 3(g. ■ n)n]d=^r - ^{ug - l)g- . (22) 

Here Ug = v{'qN) is the galactic value of v [{vg — l)g^ = (1 — Atp)gg]. Prom this follows the 
expression for the dimensionless, surrogate, anomalous quadrupole: 

'^^'^^ ^~hj l^'^^^^""' ^ ■ ""^^ ~ ^nl)\(fu, (23) 
with gjv from eq.(21), now expressed in units of Cq, and u in units of Rm- We can write in spherical 



from a potential, is not equal to g* — gw, even though the two fields have the same divergence everywhere; g* is not, 
generally, derivable from a potential. The two differ by a curl field. 
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coordinates 

= -3/ ^ dau-i)[n^p^{0 + u-'P2m, (24) 

Jo J-i 

where Pg is the is Legendre polynomial, and with = z^[(?7^ + u~'^ + 2?7^«~^^)^/^] (^ = cos^). 
Because the angular integral vanishes when the factor — 1 in the integrand is ^ independent, we 
can replace the 1 in i/ — 1 with a function of u alone, without changing the result. For example, we 
can replace this 1 with the value of v for ^ = 0, i.e., with z^[(ry^ + u~'^)^^'^], or by the vahic of v at 
^ = ±1: z^(|7?jv i M^^l). With these, convergence is improved also at large u, where v — >■ z^(?7Ar) 7^ 1. 
Changing variables to = it~^ we can also write 

/•OO rl 

q{ri) = -3 y^ dv j^^ di{v - + v^P2m, (25) 

with V = v[{r]l + v^ + 2r]i,v'^^Y/\ 

In Appendix C, I consider further one particular case: the limiting form /x(x) = x for x < 1, 
and ;u(x) = 1 for X > 1. This is useful pedagogically, since it lands itself to further approximation 
that will lead to a closed expression for the surrogate anomaly. This can be used, for example, as 
a touchstone for checking numerical results. It also affords a simple description of the phantom 
matter; and, it also seems to give a lower limit to the anomaly for a large class of interpolating 
functions. 

For this limiting form of with r^jy > 1 (which is relevant for the SS), we have a well demarcated 
region V outside which pp = 0: We have at large distances in all directions |V^| |gg| > ag, which 
takes us outside V; so V is bounded in all directions. The boundaries of V can then be obtained 
in closed form, and are described in B.7. For 77^ < 1, or for choices of ijl{x) that reach the value 1 
only asymptotically at high x, the region V is not bounded. 

The further approximation that I apply to the limiting case can be justified in the limit of 
large 1] = rj^, and leads to [see eq.(C8)] 

q(r]) = -r/-'^/V30. (26) 



^^This expression can be used when 6i^(x) vanishes faster than x~^^^ for a; —> 00 (equivalent to i' — 1 vanishing 
faster than Otherwise, the u integral diverges and we are not justified in choosing this order of integration. 

For example, for the isolated, spherical case (?7jv = 0), this expression vanishes for any choice of //, which is a wrong 
result when (5/i(x) vanishes as x~^^^, or slower. 
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(28) 



5. Results for the dimensionless quadrupole anomaly 

I use eq.(24), or eq.(25), to compute the dimensionless, surrogate quadrupole anomaly q{r]) for 
various forms of /n(x) and different r] values.^^ I also wrote a simple nonlinear-Poisson solver (based 
on the detailed blueprints described in Milgrom 1986a) to calculate q for several control cases from 
the exact expression (14). 

I consider three one-parameter families of interpolating functions, defined and discussed in 
Milgrom & Sanders (2008), and checked to give good results in rotation curve analysis in the 
transition region. They also have a rapid decline of 5ijl{x) at high x, with a negligible anomaly 
resulting from the finiteness of 5^{x). The first family is 

/^«(-) = (TtI^' 

with the corresponding 

[ l + (l + 4t/-")V2 

My) = 2 

which I use with a > 3/2. The limiting form of fj,{x) corresponds to very large a. 
The second and third families are defined by their z/(y) functions: 

I>«(y) = (l-e-^)-V2 + «e-^ (29) 

and 

u^{y) ^ (1 - e-^")(-V2«) + (1 _ l/2a)e-^", (30) 

with the corresponding /ia(x), and fla- The choice of a permits a very fast transition to the New- 
tonian regime near y = 1. These choices are only indicative; they are neither fully representative, 

nor exhaustive. 

Table 1 gives the numerical results for ~q. For some entries I also give in parentheses the value 
obtained from the numerical solution of the modified Poisson equation via eq.(14). 

The value of rj for the galactic field in the solar neighborhood is given by 

'^^^■^2201.1^ s-i) (sK^) (lO-8cm s'^) ' ^^^^ 



^^As a check of the numerics I also use eq.(22) to calculate the z component of the anomalous acceleration along 
the z axis, near the sun, to ascertain that is vanishes at the sun's position, and to check how good the quadrupole 
approximation is away from the sun. Both aspects check well; for example, the departure from a pure quadrupole 
behavior is, as expected, of the order or u/Rm, only about one percent within 80a.u.. 
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with Vq the orbital velocity of the Sun in the Galaxy, and Rq its orbital radius. ' The relevant 
range for our problem is thus roughly ^ i] < 2. In the table this gives a range 0.01 < —q < 0.3. 
The interpolating functions in the last two rows are rather less likely; if we omit those the range 
shrinks to 0.07 < -q < 0.3. 

We see that the surrogate approximation works rather well, and by and large produces q 
values that do not differ from q by much. We also see that that expression (26) is a rather good 
approximation, in the limiting case, down to 77 values of 1.5, even though it can be justified as an 
approximation only for large values of rj. 

Notice that \q\ decreases for increasing rj, when 77 3> 1, because the phantom matter becomes 
less and less important; \q\ also decreases with decreasing 17, for 77 < 1, because the phantom matter 
becomes more and more spherical. The maximum effect seems, in fact, to occur for the range of rj 
values around that for the galactic field near the sun. 

6. The effect of the anomalous quadrupole on planetary motions 

For |u| <^ Rm, the quadrupole anomaly due to the MOND EFE is of the same form as that 
produced by tidal effects of a static distant mass^*^ (such as a molecular clouds, a nearby star, etc.). 
For example, a solar mass star at a distance d is expected to give q*^ « —4 • 10~^((i/pc)~^Oo/-RM ~ 
—4 ■ 10~^(g/0.1)~^((i/pc)~^(/"^. The tidal effect of the Galaxy is dominated by the gradient in the 
direction perpendicular to the disc, call it the x axis, and we have for its quadrupole component 
q*J^^ 47rpo, where po ~ 0.076Mq/pc^ (Creze & al. 1997) is the midplane dynamical density of the 
Galaxy at the sun's position. We thus have q^J^" ^ 3 • 10~^(g/0.1)~^g^^. The other two quadrupole 
components of the galactic tidal field are about an order of magnitude smaller. Such effects are 
thus much smaller than the EFE quadrupole anomaly. 

One may wonder what the MOND effects are of the planets on each other. Beside the standard 
physics, mutual effects, which are always taken into account, there are new MOND effects here as 
well. Again, some of these have to do with the remaining small value of l — fi{x) at high x, which we 
ignore here. The effects that remain even with the limiting form of p can, again, be described as the 
Newtonian effects of the phantom matter. For any given configuration of the SS there are several 



^ The value of ao is determined from several of its appearances in the phenomenology to be around 10^ cm s^ . 
For example, Begeman, Broeils, & Sanders (1991) found from rotation curve analysis ao w 1.2-10~*cm s~^. Recently, 
Stark, McGaugh, & Swaters (2009) found that gas dominated galaxies satisfy the mass-asymptotic-rotational-velocity 
relation predicted by MOND-aka the baryonic Tully-Fisher relation-(M = aQ^G~^V^) with Oo = 1.2 • 10~^cm s~^, 
and a range of (0.7 — 2) • 10~*cm s~'^. 

^^When dealing with multi-field theories rather different r; values may appear. 

^^At radii of order Rm or larger the effect is very different from tidal effects, being the Newtonian field of pp. So, 
for example, it decreases beyond ~ Rm- 
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critical points, each begetting its own cloud of phantom matter; for example, one near each planet 

due to its competition with the Sun (Bekenstein & Magueijo 2006 consider the potential for probing 
such regions directly). There is also some distortion of the main phantom matter distribution near 
the solar-galactic critical point. I intend to discuss the problem of a few-body system elsewhere. 
Here I only mention that these effects are negligible in the present case. 

To test for the presence of the quadrupole anomaly, and constrain its parameters, one needs 
to incorporate the anomalous force in the equations of motion in programs that fit for the observed 
motions of the planets. Today, the direction to the galactic center, which we take to define our 
symmetry axis, is only about 6 degrees away from the ecliptic. To first approximation, it can thus 
be assumed that our symmetry axis is in the planetary orbital plane, and take the x axis to also 
lie in this plane.^^ Then, the anomalous acceleration in the orbital, x — z plane, is 

g» = |i(x,-2z), (32) 

at position u = {x, z) from the sun. 

To get an idea of the size of the expected effects of the anomaly on the planetary motions, 
consider, for example, the anomalous perihelion advance rate it produces. Let 

u = -, [66) 

1 + e cosV' 

be the unperturbed Keplerian orbit of a planet, with the azimuthal angle ip (which increases with 
time) measured from perihelion; (a and e are the semi-major axis and eccentricity, respectively). 

One can use the Gauss equations to calculate the variation of the orbital elements under the 
action of a given acceleration field (Brouwer & Clemence 1961). For this, one needs the instanta- 
neous radial, tangential, and normal components of the acceleration, g^, g^, (in the direction of 
increasing ip), and g^, respectively. In our case they are 

5r" = ^[l-3cos2(V-V'o)], (34) 

9? = ^75^sin(^ - ^o)cos(V' - V'o), (35) 

and I take to vanish (in the ecliptic it is about an order smaller than the other components). 
Here ipo is the value of ip on the positive z axis. With g^ = 0, there is no precession of the orbital 
plane, and the instantaneous perihelion advance, for example, can be calculated as 

= ^ -^[-g^cos^ + 5?(1 + -)sinV], (36) 

nea p 



^'^This is not a good assumption for Pluto and Icarus, which I also consider. It also breaks down at other phases 
in the galactic rotation. 



where lv is the longitude of periheHon of the perturbed trajectory, n = (MaG/a^)^/^, and p = 
a(l — e^). To get the mean periheHon advance rate we average this rate over a period. I find 




qao 



K(e)A(V'o) 



(37) 



nR 



■M 



where 



with s = [1 




(38) 



So, 



<.,«4.4..0-V„A„.,(g)(,^)-(^)-s- 



-1 



(39) 



Note that, depending on ipo, the effect can take up both signs, as X{tpo) can vary between —3/2 
and 9/4. 

I also calculated the predicted rate of secular changes in the semi-major axis due to the 
quadrupole anomaly and found that it vanishes for all values of e and tpQ. 

The kX values, and the predicted anomalous perihelion precession rates for the planets, and 
for Icarus, are shown in Table 2. 

Also shown, under |Au;|, are some indicative intervals for u from published best-fit results 
for the precession rates, obtained through fitting to specific conventional models. They are shown 
as an indication of the level of accuracy that has been reached for the different planets. The 
tightest values are for Earth and Mars and, for our reference value of g = 0.1, are ~ 8 — 10 times 
larger than the expected rates for the two planets. So, parameter values that give q ~ 0.3, and 
Oo = 1.2 • 10~^cm s~^, could already be only a factor 2-3 below these intervals. 

For the outer planets, the effect is expected to be larger, but so are the measurement errors, 
at present. So the prospects of first constraining the anomalous MOND EFE would seem to come 
if the measurements for the outer planets can be improved. 

It is interesting to compare the above results with what we get from the harmonic anomaly 
predicted in the inner SS-even without an external field-if the asymptotic behavior of the inter- 
polating function is fx{x) — 1 — Ax~^^^. In this case we get generically in MOND, an additional 



^®Fienga & al. (2009) quote Pitieva as giving a value of (—6 ± 2) • 10~^"/c, for the residual perihelion precession 
of Saturn, which might be viewed, formally, as a positive detections of an anomaly. Other analyses they quote give 
(—10 ± 8) • 10^''"/c. Such an anomaly could be explained, e.g., by the harmonic anomaly I discuss below. However, 
since these authors seem not to view these results as definitely indicating a detection of an anomaly, I took it as an 
indicative accuracy of |Ad;| ~ 10~^"/c for Saturn. 
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anomaly near the sun, with 

= (40) 

and with = 0. It is thus of the same order, and has the same distance dependence as the 
quadrupole anomaly. Using these in the Gauss equations we get a secular perihelion rate: 



{J-} ^ -6.5 . 10-«^«(e) (^)"%-' 

\10"^cms~^/ Va.u./ 



or. 



. lO^^^cm s^-^ J Va.u. 

I give in Table 2 the predicted rates for this effect for the value A = 2/3, which corresponds, for 
example, to the form //3/2 = + x^/^)~^/^, which has the appropriate asymptotic form. 



-4.1 . lO-'A.ie) (^f'Vc. (42) 



7. Discussion 

In the framework of the modified-Poisson formulation of MOND, one expects an anomalous 
quadrupole acceleration gf{u) = —q^ju^, at position u with respect to the sun, with qfj diagonal, 
axisymmetric, and traceless: —2q^^ = —2qyy = g"^, = q{r]){aQ / Rm) ■ Interestingly, this anomaly 
remains finite even if the MOND interpolating function is arbitrarily close to unity at high acceler- 
ations, despite the fact that the accelerations near the sun arc S> a,^. I find numerically \q\ ~ 0.1, 
and this corresponds to an anomalous acceleration ~ 10~^ao for the inner planets, but ~ 10~*^ao 
for Saturn, and even larger for the yet farther planets. These do not conflict with published SS 
bounds on anomalous perihelion precession, but useful constraints seem to be within reach. 

This effect can be conveniently envisaged as the tidal effect on the solar system, of the imaginary 
phantom mass that, according to MOND, hovers around the sun. 

In all the examples I studied, q turned out to be negative (repulsion from the Sun along the z 
axis and attraction in the perpendicular direction). I was not able to prove that this is always the 
case from basic properties of ji. This question of the sign of q requires a more systematic, extensive 
study. 

Our results here also apply to the unusual sort of gravitationally polarizable matter proposed 
by Blanchet &; Le Tiec (2008, 2009) as a model for MOND, as it reproduces the potential of the 
nonlinear Poisson equation. 

I find that the strength of the effect is sensitive to the form of the interpolating function ijl{x) in 
the MOND-Newtonian transition region (a; ~ 1). The constraints we have on this aspect of n comes 
mainly from rotation-curve analysis. Different studies-e.g. by Famaey k, Binney (2005), Zhao & 



-19- 



r] 


0.5 


1 


1.5 


2 


3 


5 


10 


U2 


3.4 • 10- 


-2 


0.12 




0.20 


9.8 • 10-2 


2.5 • 10- 


-2 


3.8 • 10- 


-3 


3.1 


• 10- 


-4 




3.5 • 10- 


-2 


0.11 




0.19 


0.16 


6.3 • 10- 


-2 


1.1 • 10- 


-2 


9.2 


• 10- 


-4 




3.6 • 10- 


-2 


0.11 




0.17 


0.21 


0.18 




9.0 • 10- 


-2 


1.2 


• 10- 


-2 








(0.16) 




(0.21) 


(0.24) 


(0.20) 














Ai 


3.2 • 10- 


-2 


0.11 




0.21 


0.28 


0.28 




0.13 




1.8 


• 10- 


-2 




3.9 • 10- 


-2 


8.8 • 10- 


-2 


0.11 


0.12 


0.12 




0.11 




8.2 


• 10- 


-2 








(0.10) 




(0.12) 


(0.12) 


(0.12) 
















4.0 • 10- 


-2 


8.1 • 10" 


-2 


7.9 ■ 10-2 


6.5 ■ 10-2 


4.4 • 10" 


-2 


2.3 • 10- 


-2 


8.8 


• 10- 


-3 


A* 00 


4.3 • 10- 


-2 


4.4 • 10" 


-2 


8.8- 10-3 


3.1 • 10-3 


7.3 • 10" 


-4 


1.2 • 10- 


-4 


1.1 


• 10- 


-5 








(5.3 • 10" 




(1.1 • 10^2) 


(4.6 • 10-3) 
















eg.(26) 


4.3 




3.3 • 10" 


'2 


8.1 • 10-3 


2.9-10-3 


7.1 • 10" 


-4 


1.2 • 10- 


-4 


1.1 


• 10- 


-5 



Table 1: Values of the dimensionless, surrogate, quadrupole anomaly —q for various choices of 
the interpolating function and various values of tj. In parentheses are shown some values of —q 
calculated using numerical solutions of the exact modified Poisson equation. 
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Table 2: Characteristics relating to the perihelion precession rates produced by the EFE quadrupole 
anomaly, and the harmonic anomaly: k\ as defined in eq.(38); the predicted precession rate {d>) 
for the quadrupole anomaly (calculated assuming that the orbit is in the ecliptic-not a very good 
approximation for Pluto and Icarus.); the mean anomalous quadrupole acceleration on the orbit, 
\g'^\ = \qzz\o-', the precession rate due to the harmonic anomaly (w^); its mean acceleration \g^\ = 
Aaott/RM', and an indication of existing tightness of conventional fits, |Au;|, that I was able to find 
in the literature: from (a) Pitjeva (2005), (b) Pitieva as quoted in Fienga Sz al. (2009), (c) Fienga 
&; al. (2009). All precession rates in units of 10-^" /c All accelerations in units of 10--'^2c]ji g-2 
All values are calculated for Uq = 10-^cm 3-2, q = —0.1, and A = 2/3. 
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Famaey (2006), Sanders & Noordermeer (2007), and Milgrom k. Sanders (2008)-have shown that 

rotation-curve analysis can, at present, be used to constrain ^ in its transition region. However, 
all existing studies make use of the prediction of modified-incrtia formulations. These establish 
a universal, algebraic relation of the form (1) between the Newtonian and MOND accelerations 
for circular orbits (Milgrom 1994). For consistency, such attempts to constrain for use in the 
present context should use the predictions of the modified Poisson formulation (or a multi-potential 
theory, if we are dealing with one). This is not an easy task since this theory does not provide an 
easy-to-use formula, but would require for each galaxy model, and for each choice of parameters, a 
separate numerical solution of the nonlinear Poisson equation. Brada & Milgrom (1995) did show 
that the differences in prediction of the two formulations are not very large, but they are of the 
same magnitude as the differences produced by different interpolation functions. 

Interesting constraints on the form of in the transition region will also come from studies 
of the dynamics perpendicular to the galactic disk at the solar position (Milgrom 1983, and see a 
recent study of the prospects of such analysis in Bienayme &; al. 2009). 

We thus cannot, at present, predict the exact strength of the effect, for lack of exact knowledge 
of the interpolating function and the value of 77. However, the presence of the effect can be easily 
identified if we can observe the anomaly in the motion of more than one planet: The effect for all 
planets hinges on only one parameter, q, and the dependence of the effect on orbital parameters is 
very distinctive. The anomalous acceleration, and the precession rate, increases with orbital radius 
in a prescribed way. The precession rate depends also on the orientation of the orbit relative to the 
galactic center in a known manner [through A(V'o)], which distinguishes it from a spherical anomaly 
(or from an aligned, axisymmetric one). In contrast, a central quadrupole, such as is produced 
by an oblate sun, for example, has very different characteristics: its effect decreases sharply with 
increasing orbital radius, and the perihelion precession rate it produces does not depend on the 
orientation of the orbital major axis with respect to the quadrupole axis (our 'ipo) when that axis 
is in the orbital plane (or perpendicular to it). 

The MOND anomalous acceleration becomes of order at radii ~ Rm and beyond^^, and 
could have important effects on object in the outer solar system, such as the long period comets. 
MOND could introduce completely new insights into the study of the formation, structure, and 
evolution of the Oort cloud, and shed new light on possible mechanisms for directing Oort cloud 
object into the inner solar system. But considerations of these sort are beyond the scope of my 
present discussion. 



At these radii, the field no more has the quadrupole form that which obtains near the sun. 
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A. Properties of the phantom mass 

There are several properties of the phantom mass that we can derive analytically. 

A.l. The New^tonian field of the phantom mass vanishes at the position of the sun 

This is a simple corollary of a more general result: If p{r) is a confined density distribution of 

total mass M, in a constant, background (MOND) field g^, then the force acting on this distribution 
is Mgg (Bekenstein and Milgrom 1984). On the other hand this force is also Mgg plus the force 
due to the phantom mass. It follows then that the total anomalous force on the true mass p always 
vanishes exactly. This is a generalization to the case of an external field of the following observation: 
A mass distribution does not self propel in either Newtonian gravity or in the modified-Poisson 
formulation. This means that the difference of the net MOND force and Newtonian force on the 
totality of (true) mass must vanish. But this difference is just the Newtonian force of the phantom 
matter on the true matter, which must then vanish. The phantom matter does, of course, exert 
differential forces on different parts of the true mass-thcsc arc the MOND corrections to the inter- 
body forces in the system-but they all add up to zero. In our present configuration the force on 
the Sun is MQgp(O), where gp(0) is the Newtonian field of the phantom mass at the position of 
the sun, and it must thus vanish. In the very near vicinity of the Sun the field of the phantom 
mass-i.e., the anomaly we are after-is predominantly a quadrupole. 



A. 2. The column density along the symmetry axis 

On dimensional grounds, the phantom column density along the z axis can be written as 

Sp(0) = ^7W- (Al) 

Integrating expression (10) for pp along the symmetry axis in three segments between which e 
changes sign: from — oo to the critical point, from there to the sun's position, and from there to 
oo, we get 

poo 

7 = U{oo) - W(0) = / L{x)dx (A2) 
Jo 

[L(0) = 1, L(oo) = 0]. The dimensionless factor 7 is thus independent of 77 and depends only on 
the choice of p. For the limiting form of p{x) we have U{oo) — U{0) = 1, hence in this case the 
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column density is exactly ao/^irG. For 112 we have 7 = 7r/2. Note also that 7 diverges if 5ii{x) 
behaves at large x as or slower. 

A. 3. The phantom density takes up both signs near the critical point 

At the critical point \'^(t)\ vanishes, so it must be increasing in all directions emanating thence. 
Take any surface, S, of constant |V(/>| that surrounds the critical point, but excludes the sun. 
Applying Gauss's theorem to the modified Poisson eq.(2) for the volume enclosed by S, which is 
devoid of (true) mass, gives J^, V0 • du = 0. Thus, since V|V^| is parallel to dcr everywhere on the 
surface, V|V0| • must take up both signs on every such surface. Expression (9) for pp then tells 
us that so must pp. 

A. 4. The total mass 

The total phantom mass is Mp = MQf3{ri). The dimensionless factor, /?, can be found by 
applying Gauss theorem to eq.(8) in a spherical volume whose surface S tends to infinity. 

Mp = ^jj,l-p)V<t>-da. (A3) 

Here is the MOND field; its asymptotic form is known analytically and given in eq.(32) of 
Bekenstein & Milgrom (1984). Using this expression one finds 

/ r \ 1/2 

m = M,-^L-V2sin-i I^^J - 1, (A4) 

where Lg = L{r]), pg = p{r]). 

With the limiting form of p, and 77 > 1, we have pg = 1 and Lg = 0; so the total phantom 
mass vanishes (but not pp, of course). It is easy to see this directly. In this case V is bounded, and 
excludes the sun. Outside its boundary p = 1. So, on any boundary, S surrounding V we have 
pV4) = Vcf). The surface integral of /xV0 vanishes, if E excludes the sun, as it is proportional to 
the true mass in V, and this means that Jj, • dcr = AirGMp also vanishes. 

If p does not have the limiting form, Mp is finite. For example for p = p2 = x/{l + rc^)^/^, 
we have for 77 = 2, /? 0.05, and for = 1, /3 Ri 0.23. The phantom mass is always non- vanishing 
if p{r]) < 1. But note that Mp in itself is not an indication for the strength of the anomalous 
quadrupole near the sun. Taking an extreme case, when r/ ^ 0, we have Mp 00 (reflecting the 
fact that it then constitutes an asymptotically isothermal phantom halo around the central mass). 
But the quadrupole anomaly actually goes to in this case as the phantom mass forms a spherical 
cavity around the sun. 
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AU the above applies to a single star in a constant field. In reality, the phantom masses of a 
star cuts off as it runs into those of other stars. It is easy to see that the phantom mass of a galaxy 
is just the totality of phantom masses of individual constituents; e.g., by applying Gauss theorem 
to a large volume written as the union of its sub volumes. One can then mimic the effects of MOND 
with a smoothly distributed dark matter, but only when the dynamics is probed with very course 
graining such as in rotation curve analysis. Unlike the hypothetical CDM, the MOND phantom 
matter is highly granular on the scales of distances between object in the galaxy, with regions of 
negative density, etc.. This, including the EFE we discuss in this paper, are not mimicked by CDM. 



A. 5. Space integrals weighted by pp 

One encounters various quantities of the form 

F = j pp(r)/(r)rfV, (A5) 

where /(r) can be any tensorial function of r. Examples are the various multipoles of the phantom 
mass, and its Newtonian potential and acceleration fields. Here I derive useful expressions for such 
quantities writing them as integrals over itself, instead of pp, which involves second derivatives 
of 0. These are quite useful for numerical evaluations. 

Using expression (8) for pp we can write 

P = ^J • [(1 - /^)V0]d3r. (A6) 

Integrating by parts and applying Gauss's theorem 

^ = 4^ / ^^""^^^ " ^^^"^ ■da-^j{l-p){y4>- V)f{r)d'r. (A7) 

We have already had one example of such an expression where F is the total phantom mass, with 
/ = 1. The second integral then vanishes, and the surface integral at infinity can be evaluated 
for the case of a finite external field (for g^ = the surface integral diverges). In general, we do 
not have a prescription for calculating the surface integral. However, there are cases when we can 
evaluate it. For example, if p takes its limiting form, and if 77 > 1, the surface integral vanishes 
identically over any surface surrounding V, and we can write 



In this case we get for the dipole moment, with / = r 

T> = ^J^{p-l)Vcpd'r. (A9) 
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For the quadrupole we have from eq.(A8), with / = —{r^Sij — 3rirj)/2, 

where e is a unit vector in the direction of V^. In particular 

Qzz = - l)(r • W - 3z(P,,)d'r. (All) 

The anomaly itself given in eq.(ll), corresponds to / = G|r — R|~^(r — R). In this case, 
because of /(r) ^ at infinity, we can evaluate the surface integral in eq.(A7). We have for any 
choice of //(x) compatible with MOND, and for any r] 

g«(R) = i- I ^r^[V</. - 3(V</. • n)n]d'r - 1(1 - ^,)g„ (A12) 

where n is a unit vector in the direction of r — R. 

Another version, which might also be useful, is gotten when writing expression (8) for pp as 

= ■ t"^*^ ~ W'"^)^'^] + [( V/^.) - 1]P- (A13) 

This now gives 

F = (1//X, - l)Ff> + (m/Ms - 1)(V.^ • V)f{r)dh " / (-"/^"^ " l)/(r)V^ • da, (A14) 

where is the quantity F calculated with the true mass, and the convergence of the surface 
integral has now been improved. 



A. 6. Asymptotic behavior of pp 

The asymptotic behavior of pp can be derived from the known behavior of the MOND field 
there. As derived in Bekenstein &; Milgrom (1984), the internal potential is, asymptotically (for 



MG / 2 2 



2 \ -1/2 



X +y + r-—r- ■ (A15) 



The asymptotic acceleration field is thus 

MG f 2 2 z 



and the phantom density 

-M La (2,2, \ ^ ( 2 , 2 
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M 



La 



u'-\l-^/'/u''), 



(A17) 



47r;.,(l + L,)3/2 

where x' = x, y' = y, z' = z{l + Lg)~^/^z, u' = (x'^ + y'"^ + z'^)^/^. We see that pp decreases 
asymptotically as u~^. 

The z component of the anomalous acceleration on the z axis and the x component on the x 
axis are respectively 



MG 



1 

1 



MG 



x^ 



(A18) 



We learn from this that at the same distance along the two axes, the MOND acceleration on 
the symmetry axis is larger by a factor (1 + Lg)^/"^ than that in the perpendicular direction: there 
is a stronger pull towards the Sun along the z axis. This is then also the case for the anomalous 
acceleration, since the Newtonian one is spherical. In fact, we see that the anomaly along the z 
axis is always attractive towards the sun, asymptotically. This is the opposite to what a negative 
value of q gives near the sun. 



B. Properties of the surrogate phantom mass 

B.l. Its Newtonian field vanishes at the position of the sun 

Here I show that the Newtonian field of the surrogate, phantom matter vanishes at the position 
of the sun. I first show, in general, that the force of a true density p on its surrogate, phantom 
mass vanishes identically. This force is 

= - y ppV^d^ = ^ y (V • g* + AnGp)V(p''d^r, (Bl) 

where V^'' is the Newtonian field of p alone, and the integral is over all space. The second term 
vanishes, because it constitutes the Newtonian force of p on itself. The ith component of can 
then be written as 

= 4^ / V • gVf.d^r = I {Hl),kf,d'v = -^l [(^<Ar.f.),fc - '^4>%<l>i,]d'r. (B2) 

Here, 0^ is the potential of the full Newtonian field-that of p plus the background field: gjv = 
— V^^, and v is a function of |gjv|/oo. (A subscript ,i is the derivative with respect to the ith 
coordinate.) In the second term in the integrand we can replace ^''^ ^ with 0^^^, because (j)''- and 
(1)1 differ by the constant background acceleration. We then write this term, now in the form 
'^^'k^'ki ~ ^Ai'^^'^ /^o)^], where N'{z) = z/(y^)/2. Applying Gauss theorem we then have 

F^' = --^ I {vV(t>''V(t>'' ■ da - Nda), (B3) 
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where the integration is over the surface at infinity. In the first term V^'' decreases as r~^; so 
the other factors can be taken at their constant hmiting value at infinity: v becomes Ug = v{rjj^). 
In the second term, we take the first order expansion A^[(V0^/ao)^] N{j]j^) — v^<^f ■ g^. The 
constant, hmiting value of N gives a vanishing integral, and we thus have 

= ^ / (V</.''g; • da - \I<^ ■ g^da). (B4) 

Taking the z axis in the direction of — g^, we have 

= / {^<^'<^^^ - '/'".da). (B5) 

By applying Gauss's theorem in reverse, writing the right hand side as a volume integral, we see 
that the two terms cancel, and thus = 0. This also means that the net force of the phantom 
matter on the true matter always vanishes. For our special case, it follows that the Newtonian force 
of pp on the Sun vanishes. However, since the Sun is a point mass, it follows that the Newtonian 
acceleration field of pp vanishes at the position of the sun. 



B.2. Its column density along the symmetry axis 

It is easy to gather that pp can be gotten from equations exactly like eqs.(9)(10), with V(jf) 
there replaced by — g*, with the same lA function. Its column density along the symmetry axis is 
thus identical to that of pp, and equals aoG~^7 with 7 = (27r)~^ Jq°° L{x). 



B.3. it takes up both signs 

As in the case of pp, pp oc V|g]v| • gjv, and for similar reasons this takes up both signs on any 
surface of constant |gjv| surrounding the critical point, and excluding true matter. This can also be 
seen explicitly in expression (C3) for pp, which is a good approximation for pp around the critical 
point. 



B.4. The total mass 

Use eq.(20) for pp to write 



Pp = -4^V-[(z.-l)g,]. (B6) 



Applying Gauss's theorem to eq.(20) we have 

1 

"4^ 



^v = -^ /(z.-l)g^-da, (B7) 
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on the surface at infinity, where we can write En ^ gg + s'^- Here is the Newtonian external 
field, and g'' = —MGuu~^ is the asymptotic Newtonian field of the true mass alone. Expanding 
around g^, we get 

Mp/M = v{ri^) - 1 + ri^v'{ri^)/2,, (B8) 

where r/jv = fJ'{'n)v is the Newtonian galactic field. Using the relations I'irjN) = and i^'(?7jv) = 

—rj~^ IJ,{ri)^'^ L{ri)[l + L{ri)]~^ we have 

M /M = ^ + ^ - 1. (B9) 

Since < < 1 this gives very similar values to those of the exact expression (A4) . In the MOND 
regime Lg ^ 1, f^g <^ 1 it gives compared with 7r/4/x for Mp. 



B.5. integrals weighted by pp 

We shall also need similar expression to eq.(A7) when pp is used instead of pp. Repeating the 
arguments in A. 5 using expression (B6), wc see that we get an expression similar to eq.(A7) with 
V0 replaced by — g*. For example, the surrogate anomaly is 

r(R) = ^J - 3(g^ • n)n]dh -^{ug- l)g- (BIO) 

The surrogate quadrupole moment for the limiting form of p., which I shall need below, is 

Qzz = y"(i. - l)(r • g^ - 3zg^)d'r. (Bll) 



B.6. The asymptotic behavior of pp 

Using the algebraic relation (1) it is readily seen that the asymptotic behavior of g* is 

MG ( z 



The asymptotic behavior of Pp derived from eq.(B12) is 

M L 



^*-^9^ -—3 ^' 2/' TTT- ■ 



\u = (x^ 



+ + z^)^/^]. We see a similar normalization and radial dependence as in pp, but a 
different angular dependence. 
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B.7. Delineation of V for the limiting /x and r/^v > 1 

The region V is defined by |g*| < Oq, or equivalently by Ig^rl < a^. The Newtonian acceleration 
at position u relative to the Sun is given in eq.(21). The ray from the Sun making an angle 9 with 
the z axis cuts V at 

_i_ _ r]NCos9 ± (1 - 7?^sin^0)^/^ 
~ [ (r?2 - l)cos-2e 

which give the boundaries of V- The viewing angle of V from the Sun occurs where u'^ = , for 
which sin^? = l/r]N- V intersects the z axis at = [rjf^ ^ The critical point is at rj,^^^'^ (all 

in units of Rm)- Thus, the linear size of V decreases as r/^ for large r]}^, while the distance of V 

— 1/2 

from the Sun decreases only as rjj^ . This will become important below. 

For the special = 1, V extends everywhere to the positive z side of the surface 

z_{r) = {[{r'^ + 2)^/2 _ y,2j^2}V2^ ^ith the critical point at ^ = 1. 



(B14) 



C. Further approximations for the limiting form of // 

Here I derive an approximate, closed expression for the dimcnsionless, surrogate anomaly q{ri) 
for the limiting form of It can be used to check numerical evaluations of the anomaly, and it is 
also a closed expression for some sort of a lower limit on the anomaly for > 1. The approximation 
is based on the observation that the linear size of V becomes increasingly small relative to its 
distance from the sun, as rj^f increases (see B.7). In this limit we can use for gjy in V its linear 
expansion around the critical point, where g^v vanishes. This means writing 

3/2 

gjv ~ - ^"l "° (r - SzBz), (CI) 

-TLm 

instead of the exact expression (21) (r is now measured from the critical point, unlike u, which is 
measured from the sun). With this approximation, the region V becomes the interior, P, of the 
ellipsoid 

+ + 4z^ = f, (C2) 
centered at the critical point, where i = tIj^^^^Rm- 

The phantom density in this approximation, pp, is given by 

__1 3/4 x' + y'-8z' 

Pp- qPoVn (^2 +y2 + 4^2)5/4' i^^^ 

where = SMq/AttRIj; it is axisymmetric, but also symmetric under reflection in z. This means 
that, unlike pp, and pp, it does not have a dipole moment. One can check that the Newtonian field 
of Pp does not vanish on the z axis, except at the origin. In particular, in itself, it does not vanish 
at the sun's position, as required. We cannot then use expression (CI) directly in eq.(23) to get 



-29- 



the desired approximate expression. In replacing pp by pp, the quadrupole moment of the phantom 
mass, Qij (not to be confused with the quadrupole anomaly), is expected to be well approximated, 
but we have completely lost the important contribution of its dipolc moment D, which is needed to 
balance the quadrupole and annihilate the field at the sun's position. The dipole moment derives 
from the small z asymmetry that does not affect Qij much. I thus use the above approximation 
only for calculating the quadrupole moment. I imagine that pp is slightly distorted so as to attain 
a dipole moment. If we know the quadrupole moment wc can determine the dipole moment by 
requiring that they cancel at the sun, use this to eliminate the dipole, and express the whole effect 
in terms of the quadrupole. To recapitulate, I approximate the anomaly-the Newtonian field of the 
phantom mass-by its dipole-plus-quadrupole contribution, writing for the potentiaP'^ 



l'{R) = -^T>-R--^R'R'Q,j. (C4) 



(R, like r, is measured from the critical point.) Now, to have the acceleration field of (p"" vanish at 
the sun's position we have to have^-*^ 

D = ^Qzzez- (C5) 

2ilo 

~ ~ —1/2 

(Here, Rq is the distance from the Sun to the critical point of V; Rq = rj^ Rm-) I put this value of 
D back into the expression for to get the desired anomaly in this approximation, which requires 
only knowledge of the quadrupole moment. From this we can identify the anomalous quadrupole 
and find 

& = = = ^Qz. ^ -^^iv). (C6) 

Using expression (CI) for g^v in eq.(Bll)-the formula for the quadrupole moment-then gives 
straightforwardly 

The scaling with this particular power of 77^ is understood as follows: For large rj^, the linear 
size, d, of V is (X Rm'Hn^^'^ i^i directions . The characteristic surface column density is always 
S a Oo/G; the quadrupole moment, which scales as Ed"^ must then be Q oc aoG~^ R'^rj'^^ . This 
gives finally 

q{v) = (C8) 

Note that our expressions (C7)(C8), and the approximation behind them, must break down 
for values of r] slightly above 1, or lower. Wc know, for example, that the anomaly should vanish 
in the limit r]^ 0, which relation (C8) fails completely to account for. 



^''D = Jpp(r)rd^r, and Qij = — i J pp{r){r'^Sij — 3rirj)d^r. Qij is a traceless matrix, and with the symmetry of 
the problem it is diagonal, with Qxx = Qyy = —Qzz/2. 

^^This relation does not need to hold exactly since higher multipoles may also contribute to the field value at the 
sun, but it should hold to the same accuracy as the dipole-plus-quadrupole approximation. 
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